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Abstract. We present a high-resolution (R « 16,000) spectrum and a narrow-band image centered on the [Nen]12.8 /jm line 
of the central kpc region of the starburst/Seyfert 2 galaxy NGC 7582. The galaxy has a rotating circum-nuclear starburst disk, 
shown at great detail at a diffraction-limited resolution of 074 (« 40 pc). The high spatial resolution allows us to probe the 
dynamics of the [Ne n] gas in the nuclear regions, and to estimate the mass of the central black hole. We construct models of 
gas disks rotating in the combined gravitational potential from the stellar bulge and a central black hole, and derive a black 
hole mass of 5.5 x 10 7 M G with a 95% confidence interval of [3.6,8. 1]xl0 7 M Q . The black hole mass combined with stellar 
velocity dispersion measurements from the literature shows that the galaxy is consistent with the local Mbh-o". relation. This 
is the first time that a black hole mass in a galaxy except our own Milky Way system has been estimated from gas dynamics in 
the mid-infrared. We show that spatially resolved mid-infrared spectroscopy may be competitive with similar techniques in the 
optical and near-infrared, and may prove to be important for estimating black hole masses in galaxies with strong nuclear dust 
obscuration. The high spectral resolution allows us to determine the heliocentric systemic velocity of the galaxy to between 1614 
and 1634 kms~'. The mid-infrared image reveals several dense knots of dust-embedded star formation in the circum-nuclear 
disk, and we briefly discuss its morphology. 
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1. Introduction 

Much effort is currently being put into understanding how 
black hol es form and grow in the context of galaxy forma- 
tion (e.g. lHaehnelt & Kauffman dbOOOt iMarconi et alJ [2004; 
iDi Matteo et alJl2005t ICattaneo et alJl2005l). The M^7 -cr t re- 
lation dFerrarese & Merritt 2000; Gebhardt et al1 l2000l) . being 
a tight correlation between black hole mass and stellar veloc- 
ity dispersion, is strong evidence for a causal link between the 
growth of the black hole and the galaxy bulge. The Mbh-c* re- 
lation has become our most important tool for studying the co- 
evolution of black holes and galaxies. To advance our knowl- 
edge we need to investigate how black hole mass and bulge 
velocity dispersion relate in a variety of galaxies, over a wide 
range in Hubble type, black hole mass and velocity dispersion. 
However, direct measurements of black hole mass in galaxies 
are challenging because resolutions of the order of the radius 
of the sphere of influence of the black hole are required, a scale 
which is typically «1". 

A widely used method for black hole mass measure- 
ment is that of spatially resolved spectroscopy of rotating 



gas disks in the centers of nearby galaxies. Because spa- 
tial resolutions of typically 071 are needed, this method 
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2003tlAtkinson et al.l2005l) . With the superior resolution of the 
HST, gas rotation close to, or inside, the black hole's sphere 
of influence, i.e. inside which the black hole potential domi- 
nates over the bulge potential, can be measured. For the more 
massive black holes, say > 10 s M , or for sufficiently nearby 
galaxies, the sphere of influence can also b e resolved from the 
ground at near-infrare d wavelengths (e.g. Marc oni et al J200lt 
Tadh unter et all2003l) . 



Now that the Space Telescope Imaging Spectrograph 
(STIS) on board HST is no longer operational, it becomes 
necessary to find alternative techniques with sufficient spa- 
tial resolution to probe dynamics of gas and stars close to 
the black hole's sphere of influence. Adaptive optics assisted 
spectroscopy in the near-infrared is the most obvious alterna- 
tive because resolutions of < 071 can routinely be reached 
dHaring-Neumaver et alJ2006T) . Here we show that also diffrac- 
tion limited spectroscopic imaging in the mid-infrared at 8m 
class telescopes is a viable technique. 
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Some galaxies have nuclei that are hidden behind large 
amounts of dust, hence the nuclear regions become unavailable 
at optical, and sometimes also at near-infrared, wavelengths. 
However, the dust is penetrated at mid-infrared wavelengths, 
and mid-infrared spectroscopy at high spatial resolution can 
therefore be used to constrain black hole masses in dust en- 
shrouded galaxy nuclei. The starburst/Seyfert 2 galaxy NGC 
7582 is an example of a such a galaxy. Here we present data of 
this galaxy taken with the ESO mid-infrared instrument VISIR 
(VLT Imager and Spectrometer for mid-InfraRed). The data 
consist of a narrow-band image and two high-resolution long- 
slit spectra centered on the [Nen]12.8/im line. 

NGC 7582 is a highly inclined barred spiral galaxy 
(SBab) at a systemic velocity cz ~ 1620 kms~' with 
a composite starburst/Seyfert 2 nucleus. The nucleus has 
a we ll-define d one-sided ionizati on cone in [Oin]5007 
JStorchi-Bergmann & Bonattolll99l() . In the center, there is 
circum-nuclear Ha emission distributed in a rotating kpc-scale 
disk in the plane of th e galaxy, perpendicular to the cone 
axis JMorris et al.lll985 | ). The nucleus is h idden behind large 
amounts of dust (Regan & Mulchaev 1999), presumably asso- 
ciated with the putative torus. 

The starburst in the center is associated with the nu- 
clear kpc-scale d isk see n in low-exci t ation emission 

lines JMorris et all Il985t ISosa-Brito etaD l200ll) and in 
the soft X-rays, extend e d em ission from the starburst 
is seen jLevenson et alJ l200ll) . Low-resolution spectra 
and images in the mid-infrared show PAH (Polycyclic 
Aromati c Hydrocarbon) emissi on from star forming re- 
gions ( Siebenmor gen et all 12004). The nuclear AGN (Active 
Galactic Nucleus) point source is strong and well-defined 
in the L (3.8//m) and M (4.7jum) bands, and the star form- 
ing kpc-sca l e disk also shows up in the L-band image 
JPrieto et all |2002), probably because of PAH emission at 
3.3/u m. Near-infrared spec tr a show both broad an d narrow 
Bry jReunanen et a l .1 120031) . IWold & Gallium! ({2006) discuss 
the star formation in the disk in more detail, in particular the 
presence of embedded star clusters discovered in the [Nen] 
narrow-band image, with no counterparts detected at optical or 
near-infrared wavelengths. They estimate that the two brightest 
sources in the disk may contain roughly 1500 O stars each, 
and are likely to be young, super star clusters embedded in 
their own dusty birth material. 

Throughout the paper, a Hubble constant of Ho = 70 
kms~' Mpc 1 is used, hence one arcsec corresponds to « 108 
pc in the rest-frame of the galaxy. 

2. Observations and data reduction 

A narrow-band image and two long-slit spectra were obtained 
with the VISIR instrument mounted at the Cassegrain focus 
of the UT3 Melipal telescope. The VISIR detector is a SiAs 
256x256 DRS array, and the pixel scale used during the obser- 
vations was 0'.'127, yielding a field of view of 32'.'5x32'.'5. 

The image was taken on 2004 September 29 through 
the [Nen] narrow-band filter centered on 12.8 1/im (half-band 
width 0.21/im). A perpendicular chopping-nodding technique 
was used with a chop throw of 12'.'0 and a total integration time 



of one hour. The two high-resolution spectra were obtained on 
October 2 using the [Ne n] 12.8/im long-slit echelle mode. Both 
spectra are centered on the nucleus. One spectrum was obtained 
with the slit approximately parallel to the major axis of the 
circum-nuclear disk and the other with the slit perpendicular 
to the major axis. A parallel chop-nod technique was used with 
nodding along the slit. The chop throw amplitude was 12'.'0 and 
the total integration time for each spectrum was one hour. 

The spectra were obtained with the 0V75 slit and the grat- 
ing tilted to a central wavelength of 12.8803/mi. With this set- 
ting, the dispersion is 1.372xl0~ 4 /v m P er pixel, and the resolv- 
ing power R « 16000. For calibrations, a telluric standard star, 
HD1522, with spectral type K1.5III, was observed. Conditions 
on both nights were clear with humidity below 10%. 

The PSF in the narrow-band image has a FWHM (Full 
Width at Half Maximum) of (X'4. The spatial resolution in the 
spectrum is also 0V4, hence the data are diffraction limited as 
the diffraction limit for the VLT is / /39 at 12.8/mi. The spec- 
tral resolution is given by the width of the function which is 
a convolution of the Line Spread Function (LSF) and the slit 
function. If we assume that the slit function is a rectangular 
function with width equal to the slit width of 5.9 pixels, and 
convolve it with a Gaussian LSF with a FWHM of 4 pixels, we 
obtain a profile with a width of ~ 7 pixels, corresponding to 
the FWHM of sky lines in the spectra. The LSF can therefore 
be approximated by a Gaussian with a FWHM of 4 pixels, im- 
plying that the instrumental broadening (FWHM) is 4 pixels, 
or 13 km/s. 

Pipeline-reduced images and spectra are used throughout 
the paper. For each telescope nodding position the average of 
all chopped frames is built and the different nodding frames are 
combined to derive the final image. The wavelength calibra- 
tion is performed by cross-correlating the position of skylines 
in the observed spect rum with that of a HITRAN spectrum 
(Rothma n et alJ Eo03l) . The pipeline wavelength calibration 
was checked by identifying the three strongest skylines with 
the corresponding skylines in high resolution {R » 130,000) 
iV-band solar spectra from the McMath Fourier Transform 
Spectrometer ( Walla ce et alJ 11994 . and an agreement to within 
less than two pixels was found. 

3. The circum-nuclear star forming disk 

In order to study the structure of the circum-nuclear disk, we 
formed a color image of the central ~ 10" x 10" by combining 
the VISIR narrow-band image with archival HST WFPC2 and 
NICMOS images taken through the F606W and F160W filters 
(programmes 8597 and 7330, respectively). The NICMOS im- 
age was obtained with the NIC2 detector, and for the WFPC2 
image, we used only the PC chip covering the nucleus. The 
colour image is shown in Fig. ^ where the WFPC2 image has 
been coded in blue, the NICMOS image (//-band) in green and 
the [Ne n] narrow-band ima ge in red. The imag e agrees well 
with the model favoured by Mor ris et"afl dl985l) based on ob- 
servations of [O m] and Ha. They suggest a rotating, star form- 
ing kpc-scale disk surrounding the nucleus with a one-sided 
outflowing ionization cone with the axis perpendicular to the 
disk plane. 
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Fig. 1. Color composite of the central kpc. The optical WFPC2 
F6Q6W image is coded as blue and the near-infrared NICMOS 
F160W image as green. The narrow -band image centered on 
the [Nen]12.8/im line is coded in red. North is up and East is 
to the left. 

There is a striking anti-correlation between the three col- 
ors in the image. The blue-coded optical emission dominated 
by the stellar continuum is obscured by a dust band stretch- 
ing across the nucleus from the South-East towards the North- 
west. Since there is some contribution from Ha to the optical 
flux the cone-structure can be seen in blue at PA w 250° open- 
ing up towards the West. The green-coded near-infrared emis- 
sion coming mainly from stars in the circum-nuclear disk is 
seen through the dust where the optical emission is obscured. 
But even the near-infrared emission is patchy because of dust 
obscuration, and in regions where both the near-infrared and 
optical emission has been extinguished by dust, the red-coded 
mid-infrared emission peeks through. 

[Nen] is collisionally excited and the 12.8/im emission 
traces gas that has been ionized by UV radiation from newly 
formed, massive stars. Several knots of emission are seen, 
associated with dense, embedded regions of star formation. 
IWold & Galliano! d2006l) discuss the possibility that the bright 
[Nen] knots may be young star clusters consisting of ~ 10 3 
stars, embedded in dusty material. The two most prominent 
knots are seen to the South of the nucleus. The rightmost knot 
is seen i n the deconvolved 11.9^m image (tracing PAH emis- 
sion) by Sieb enmorgen et all ll2004t) . but the one to the left is 
undetect ed at 1 1 .9um. I t is also very weak, or undetected, in H2 
dSosa-Brito et all200ll) . 

Using isophotal ellipse fitting, we measure the elliptic- 
ity 1 and position angle of the [Nen] disk after masking out 

1 Here defined as e = 1 — - with a and b being major and minor 
axis, respectively 



the knots. The best measurements at a semi-major axis of « 
2'.' 5 give an ellipticity of 0.47 + 0.11 and a position angle of 
158° + 4°. Assuming a thin, circular disk, the inclination angle 
is given by cos i = -, implying that i = 58°. The circum-nuclear 
gas disk therefore has, within the errors, the sa me inclination 
angle as the galaxy, 60° ±10° jMorris et all 19 85). and is con- 
sistent with lying in the plane of the galaxy. The rotation of the 
disk follows the direction of rotation of the galaxy. 

4. The velocity field in the disk 

The two [Ne 11] longslit spectra show the velocity field across 
the disk along two directions almost perpendicular to each 
other. Fig. |2 shows the full two-dimensional spectra together 
with the narrow-band image with slit positions indicated. The 
parallel spectrum (shown to the left) was taken with the slit 
positioned over the nucleus and one of the compact knots of 
[Ne 11] emission to the South. The position angle is PA = 172°5, 
which is close to the position angle of the major axis of the disk. 
The perpendicular spectrum, shown in the right-hand panel, 
was obtained at PA = 90°, i.e. along the axis of the ioniza- 
tion cone. Both spectra show continuum emission from hot dust 
in the centre. The parallel spectrum clearly exhibits rotation, 
whereas the velocity across the minor axis as seen in the per- 
pendicular spectrum is roughly constant. 

In order to derive velocity curves, apertures two pixels wide 
were extracted from the spectra. A Gaussian was fitted to each 
extracted spectrum so that the line centroid, FWHM and flux 
as a function of position along the slit could be determined. 
The results are shown in Fig. [5] where wavelength has been 
converted to radial velocity and FWHM to velocity dispersion 
(FWHM/2.35). 

The parallel chop-nod technique produces one positive and 
two negative images on the detector. The positive image con- 
tains half of the total flux and each of the two negative images, 
a quarter of the total flux. As we were concerned with deter- 
mining the position of the wavelength centroid as accurately 
as possible we chose not to combine the three spectral images, 
as imperfect registration could have introduced centroid shifts. 
Averaging the three images also does not increase the signal- 
to-noise by more than a factor of 2/ V3 (because each image 
contains only part of the total flux). 

The data shown in Fig. [3] were extracted from the positive 
spectrum which has the highest signal-to-noise. For estimat- 
ing uncertainties in wavelength centroid, FWHM and line flux, 
we utilized both the positive and the two negative spectra. The 
uncertainty in each quantity x, where x may be line centroid, 
FWHM or flux, was calculated as (|x pos - x neg \)/2, where x neg is 
the average of the measurements from the two negative images 
and Xpo S is measured from the positive image. We find an un- 
certainty in the wavelength centroid of typically 5-10 kms -1 , 
depending somewhat on the position along the slit. There is an 
additional uncertainty due to the fact that the spatial resolution 
in the spectrum is smaller than the slit width. However, this is a 
systematic offset which does not vary with wavelength, hence 
was not taken into account here. The uncertainty in the velocity 
dispersion is 3-7 km s _1 and the uncertainty in the flux, from a 
few up to a; 10 per cent. 
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Fig. 2. The two high-resolution spectra of the [Ne n] line with 
the slit positions indicated on the narrow-band image. The par- 
allel spectrum is shown to the left, and the perpendicular spec- 
trum to the right. The two regions of increased noise at * 1350 
and 1850 kms~' are due to strong skylines. Continuum from 
hot dust emission in the center is seen. 
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Fig. 3. Radial velocity, velocity dispersion and flux of the 
[Ne n] line as a function of position along the slit in the two 
spectra. The panels to the left show the parallel spectrum, and 
the panels to the right, the perpendicular spectrum. 



The rotation curve in the upper left-hand panel of Fig. [5] 
shows that the disk rotates with an amplitude of * 200 km s _1 . 
The rotation curve also appears relatively smooth with some 
irregularities close to the edges. In the blueshifted southern part 
of the disk (over pixels 26-38) a slight turn-over in the velocity 
is observed, probably indicating that most of the mass probed 
by the rotation curve is contained in the probed region. The 
rotation of the disk is such that the emission to the North is 
redshifted and the emission to the South blueshifted. 

The perpendicular velocity curve to the right shows a more 
constant velocity of ss 1650 km s _1 across the disk with a slight 
decrease in velocity toward the East. An approximately con- 
stant velocity profile is what we expect from a perpendicular 
spectrum aligned along the minor axis of a thin rotating disk. 

The middle row of panels in Fig. [3] shows velocity dis- 
persion as a function of position along the slit. Typically we 
measure dispersions of 40-80 kms~' with maximum around 
the center of the disk. There are several different factors con- 
tributing to the velocity dispersion. Apart from unresolved ro- 
tation, there are contributions from instrumental and thermal 
broadening. The instrumental broadening is ss 13 kms~' (see 



Section |2j, and the thermal broadening broadening (FWHM) 
for a 10 4 K disk * 20-30 kms 1 . This adds up to a FWHM 
of » 24-33 kms~', or <x ss 10-14 kms -1 . We measure <x > 40 
km s _1 so there must be a non-negligible contribution from bulk 
motion in the gas. The increase in velocity dispersion toward 
the center of the disk is seen in both the parallel and the per- 
pendicular spectrum, probably caused by spatially unresolved 
rapidly rotating gas close to the nucleus. 

The radius of the black hole sphere of influence is 
GMbh/vI where Mbh is the black hole mass and cr t is the 
central stellar velocity dispersion. The velocity dispersion of 
NGC 7582 listed in the HyperLed a database 2 is 157+20 k m/s, 
and using aperture correc tions by J0reensen et al. (1995) and 
the Mbh - cr* relation by Trem aine et alJ J2002I) . we estimate 
a black hole mass of * 4 x 1O 7 M0. Hence the radius of the 
black hole sphere of influence is S 07 1 . The spatial resolution 
in the [Ne n] spectrum is a factor of 3-4 larger than this, so 
the sphere of influence is unresolved. But even though the spa- 
tial resolution is not sufficient to probe gas rotation inside the 

2 http://leda.univ-lyonl.fr/ 
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sphere of influence, an estimate of the black hole mass can be 
obtained from the steepness of the rotation curve across the 
nucleus. In the following, we set up a model of the rotating 
disk to be able to estimate the centr al black hole mass. Our 
modeling follows largely that of Bart h et all J200ll) . and is the 
standard gas-dynamical method used for nuclear gas disks in 
galaxies (see also Mar coni et alJl200ll iTadhunter et al.lEo03l 
lAfkinson et alJl2005h . The gas disk is assumed to rotate in the 
combined potential of the galaxy bulge and the central super- 
massive black hole. By varying the black hole mass as well as 
other parameters such as disk inclination and mass-to-light ra- 
tio, different models can be compared to the data in order to 
find the best fit. We describe the modeling in more detail in the 
following section. 



5. The model velocity field 

5.1. The basic equations 

We assume that the circum-nuclear disk can be approximated 
by a flat, circular disk with constant inclination rotating in the 
gravitational potential W. At every radius R in the disk, the ro- 
tation is dictated by a combination of the gravitational potential 
from the stellar bulge and the central, supermassive black hole, 



The velocity along the line of sight is given by 



CM, 



BH 



R 



(1) 



The bulge potential is and the black hole point mass poten- 
tial is GMbh/R- The bulge potential dominates far out in the 
disk, but close to the center, within the black hole sphere of 
influence, there is a non-negligible contribution from the black 
hole. The circular velocity, v c , in the disk is obtained by differ- 
entiating the potential with respect to R: 



(2) 



The plane of the sky is defined by the coordinate system 
(x',y') where the x'-axis coincides with the major axis of the 
disk. We also define the slit coordinate system (x s ,y s ) with the 
y s -axis along the length of the slit. The transformation between 
the two coordinate systems is given by 

x = (x s - b) sin 9 + y s cos 9 

y = —(x s - b) cos 6 + y & sin 9 (3) 

(see appendix B bv lMarconi et alJfcoOll) . where 9 is the angle 
between the slit and the major axis of the disk. In our case, 
the slit is rotated by 9 — —14° 1 with respect to the major axis 
(i.e. 158°4, which is the position angle of the disk, minus the 
slit position angle of 172°5). The center of the disk, i.e. the 
dynamical center, is assumed to have coordinates x s = b and 
y s = 0, hence the ^-parameter enables us to take into account 
that the dynamical center may not be at the exact center of the 
slit. 

The inclination of the disk is i, and we have by deprojection 

that 



R 2 = x' 2 + 



y 



v = Vc sin 1 - 



R ' 



(5) 



From the equations above, we calculate the velocity along the 
line of sight in the slit coordinate system for any combination of 
bulge potential, black hole mass, disk inclination and dynami- 
cal center position b. By comparing the velocity curves derived 
from the different models to the observed one, the best-fitting 
parameters, including black hole mass, can be found. 

5.2. The bulge gravitational potential 

An important part of setting up the model is to determine the 
bulge gravitational potential, *¥*(R). In order to do this, we use 
the stellar density in the bulge as a tracer of the gravitational 
potential. The NICMOS F160W image is well suited for this 
because the near-infrared is a good tracer of stellar mass and 
because the image has high spatial resolution. Moreover, near- 
infrared is less affected by dust than the optical, a clear advan- 
tage here because the center of NGC 7582 is strongly affected 
by dust. 

In order to parameterize the surface brightness distribu- 
tion, we use the Multi-Gaussi an Expansion (MGE) method 
developed by Cappellari (2002). It consists of parameterizing 
the galaxy surface brightness with a series expansion of two- 
dimensional Gaussian functions which can be converted to a 
gravitational potential *F«. But since even the F160W-band is 
affected by dust, we found it necessary to make a correction 
before applying the MGE analysis to the image. This was done 
by mapping the extinction across the central parts of the galaxy 
by forming an F606W - F160W image. The absolute extinc- 
tion in the F160W-band was thereafter calculated by assuming 
a galactic reddening law JCardelli et alll98H) and a correctly- 
normalized dust-corrected image made by matching pixels in 
the uncorrected image to pixels in the corrected image that ap- 
peared unreddened. For simplicity, we assume a uniform intrin- 
sic colour in the bulge, but in reality a gradient may exist. 

In order to avoid introducing too many distortions in the 
surface brightness profile, the dust lane that runs across the 
galaxy was masked out before performing the MGE fit. The pa- 
rameters of the best fit is listed in Tabled where each Gaussian 
is characterized by a surface density I',, a dispersion cK , and an 
axial ratio q',. The luminosity in the last column is the to- 
tal luminosity of the Gaussian. Primed quantities are projected 
whereas unprimed quantities are intrinsic. The luminosity is 
calculated as 2nl' k q',crk, where we have assumed that the gravi- 
tational potential is axissymmetric so that <j' k — cr^. 

In order to take the shape of the NICMOS PSF into ac- 
count in t he MGE analysis, we also computed a Tiny Tim PSF 
dKrist & Hooklll997|) and p arameterized it with four circular 
Gaussians ( CaDDell aril2002h : 



PSFOR') = V G jt exp[-/?' 2 /(2(rf)]/(27rcr* 2 ), 



k=\ 



(6) 



(4) 



where cr* are the dispersions and the weights of the 
Gaussians. The parameters of the MGE fit to the NICMOS PSF 
are listed in Table |2] 
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Table 1. Parameters of the MGE fit to the galaxy bulge in the 
NICMOS F160W filter. 



k 






4k 


Lk 




Lo,F 1601V pc 


arcsec 




10 6 L Q- f i6ow 


1 


33307 


0.330 


0.716 


192.5 


2 


1041 


1.333 


0.872 


119.5 


3 


4969 


2.229 


0.556 


1018.5 


4 


343 


7.188 


0.556 


730.2 



Table 2. Parameters of the MGE fit to the NIC2 PSF in the 
F160W filter. 



k 


Gk 


tr* [pix] 


1 


0.251 


0.75 


2 


0.475 


1.78 


3 


0.146 


7.64 


4 


0.127 


19.51 



The gravitational potential corresponding to an MGE den- 
sity profile in the oblate axissymmetric case is given by 



m,(R,Z= 0) = GT 



jfc=0 



L k Q k (R,z = 0) 



(7) 



femsellem etal.ll 19941 ICappellari etai1l2002l) . where Y is the 
stellar mass-to-light ratio. The radius in the equatorial plane 
is R and z — because the circular velocity in the equato- 
rial plane is determined solely by the radial component of the 
force. We assume that the gas disk lies in the equatorial plane 
of the galaxy, hence the radial velocities we derive also lie in 
this plane. The factor Qf, is defined as 



Q k (R 



,z = 0)= f 
Jo 



1 exp(-R 2 T 2 /2of) 



dT, 



(8) 



(see Cappellari 2002) where is the intrinsic axial ratio, cal- 
culated under the assumption that the galaxy has an inclination 
of 56°2, as measured by the MGE fitting routine. 

5.3. Generation of two-dimensional synthetic spectra 

Having determined the bulge gravitational potential, we use the 
equations outlined in Section lSTl to calculate models for differ- 
ent values of black hole mass and dynamical center offset, b. 
The disk inclination and the mass-to-light ratio are held con- 
stant at i - 58° and T = 3.8 as the data do not have sufficient 
resolution to allow for a multi-parameter fit. The mass-to-light 
ratio is fit separately, prior to fitting Mbh and b, as explained in 
Sectionl5~4l 

Each model is "observed" with the 0775 slit and the instru- 
mental setup that was used during the observations. We also 
convolve with the resolution of the observations, and fold the 
[Ne n] narrow-band image into the model in order to produce a 
two-dimensional synthetic spectrum. We follow the procedure 
of iBarth et alJ 11200 II) . which is briefly described below. A de- 
tailed explanation can be found in Barth et al.'s paper. 



First we set up a coordinate grid on which to compute the 
synthetic spectrum. The grid is subsampled by a factor of two 
compared to the data, so the pixel size in the model spectrum 
is half the size of the real pixels. The model is therefore re- 
sampled before comparing to the data. The construction of the 
synthetic spectrum consists of calculating the observed surface 
brightness of the disk through the 0775 slit. This is done by as- 
suming that the intrinsic line profiles are Gaussians and adding 
them up at every point in the disk, weighing by the [Ne n] sur- 
face brightness (as provided by the narrow-band image) and 
smearing by the seeing PSF. 

The synthetic line profile, fy s (v), is given by Barth et al.'s 
(2001) eq. 4, and is reproduced here in order to clarify the dis- 
cussion: 



./v. (v) = J ^.V, ; /'</../\u..vJ 

[v - v p (i, j) - (x s - jc c )v d ] 



*s=*0 l,J 



x exp 



-2oi 



(9) 



The [Ne n] surface brightness in the disk at pixel (i, j) on the 
model grid is denoted S y, and the value at (x s ,y s ) of a PSF 
centered at pixel is denoted P(i, j\x s ,y s ). Hence the mul- 
tiplication of S and P represents a convolution of the surface 
brightness distribution with the PSF and takes into account 
that different parts of the disk contribute to the velocity pro- 
file at a given position. We assume that the PSF is a Gaussian 
with FWHM equal to the seeing. During the computations the 
PSF was truncated at 2.5xFWHM. Since we are effectively 
convolving the narrow-band image with the PSF, the narrow- 
band image first had to be deconvolved. This was achieved by 
block replicating the original image by a factor of two and de- 
convolving using th e Multiwavelength Deconvolution Method 
dStarck & Murta g hl2002h. 

Since the intrinsic line profiles are assumed to be 
Gaussians, the exponential factor in Eq. [9] contains the model 
velocity field v p (/, J) and a correction factor (x s - x c )vd taking 
into account the uncertainty in velocity caused by the fact that 
the recorded wavelength of a photon depends on where it en- 
ters the slit along the x s -axis. The pixel size in velocity units is 
v,j = 3.2 km s _1 , and the center of the slit is denoted x c . The ve- 
locity dispersion of the gas, <x r , occurs in the denominator. As 
seen from Fig- EJ it varies significantly across the disk, hence 
is difficult to fit by a model. We therefore approximate <x r with 
» 55 kms~' which is the average across the disk. 

In order to avoid including continuum emission in the syn- 
thetic spectra, the narrow-band image was continuum sub- 
tracted before incorporating it into the model. Because it was 
obtained solely for the purpose of aiding with the determina- 
tion of slit positions, images in the continuum filters bracketing 
the [Ne n] line were not taken. We therefore subtract a contin- 
uum based on that measured in the spectrum instead, achieved 
by first scaling a point source to the continuum flux measured 
in the spectrum, and thereafter subtracting it from the image. 
Because the narrow-band image is folded into the model, the 
position of the slit relative to the narrow-band image had to be 
determined. This was done by first rotating the image so that 
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Data Model 

12.91 [ " ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' ] [ ' I ' ' ' I ' ' ' I ' " 



12.90 - 




12.86 f i i i I i i i I i i i I i i i I i i i I i i i 1 f i i i I i i i I i i i I i i i I i i i I i i i 1 
-6-4 -2 2 4 6 -6 -4 -2 2 4 6 
Spatial axis [arcsec] Spatial axis [arcsec] 

Fig. 4. Comparison between the real spectrum (left) and a syn- 
thetic spectrum (right), calculated for a black hole mass of 
5.5 x 10 7 M and with the dynamic center shifted 0/3 from 
the center of the slit. 



the length of the slit was parallel to the y-axis of the rotated im- 
age, and thereafter collapsing the spectrum in the wavelength 
direction and matching it to the section of the image covered 
by the slit. 

The summation over x s in Eq. [9] runs across the width of 
the slit, where the edges of the slit are given by x s = xq and 
x s = x\. The y s -coordinate runs along the spatial direction 
of the spectrum. To create the full two-dimensional spectrum, 
fy s (v) has to be evaluated for every pixel along the spatial di- 
rection. An example of a synthetic spectrum is shown in Fig.|4] 
The two regions of emission on either side of the center are 
well reproduced, as are their relative compactness and bright- 
ness. As the synthetic spectrum is noise free, and also void of 
sky lines, it differs slightly from the observed one, especially 
around 12.863yum where there is increased noise due to the sub- 
traction of a sky line. Also, since the synthetic spectrum was 
produced by convolving with a continuum-subtracted narrow- 
band image, the two differ in flux at the center. 

In order to derive model rotation curves, we made extrac- 
tions from the synthetic two-dimensional spectra in a manner 
identical to that which was done for the real spectrum. 



5.4. Constraints on the black hole mass 

Given that the VISIR spectrum has a spatial resolution of « 
074, we merely fit two parameters, namely the black hole mass 
and the location of the dynamical center within the slit. The 
inclination of the disk is fixed at i = 58° and the mass-to-light 
ratio at T = 3.80. We calculate models on a grid of black hole 
mass and dynamical center location, b. The spacing in Mbh 
was varied between 0.5 xlO 7 and lxlO 7 M , and the spacing in 
b, between O.'l and O/05, depending on how close the model 
was to the global minimum in^ 2 . 

For every pair of Mbh and b a synthetic two-dimensional 
spectrum was constructed, and one-dimensional spectra ex- 
tracted from it using apertures of the same width as for the real 
data. Thereafter a Gaussian fit was made to each in order to de- 
termine the line centroid so that a radial velocity curve could be 
derived. The best-fitting model was found in a standard manner 
by minimizing^ 2 between the observed and the model rotation 
curve. Fig.|5]shows the observed rotation curve together with a 
model rotation curve extracted from a synthetic spectrum. 

In order to convert the MGE surface density profiles to a 
gravitational potential, the stellar mass-to-light ratio, T, in the 
bulge has to be known, see Eq. Q Outside the black hole's 
sphere of influence, the rotation curve does not depend on ei- 
ther the black hole mass or the exact location of the dynamical 
center in the slit. Hence the outer parts of the rotation curve are 
not sensitive to Mbh and b, but depend instead on the mass- 
to-light ratio. We therefore fit T separately, before fitting Mbh 
and b, by neglecting the inner few points on the rotation curve 
and fitting only the outer parts. Synthetic spectra were gener- 
ated with Mbh = and b — and with mass-to-light ratios in 
the range 2.4^1.8 in steps of 0.2, while keeping the disk incli- 
nation constant at i = 58°. Apertures were extracted and model 
rotation curves produced. The best-fit mass-to-light ratio was 
thereafter found by fitting the model rotation curves to the ob- 
served one with the three inner apertures (centered at pixels 36, 
38 and 40) excluded. The best-fit model has T = 3.80 M /L 
in the F160W filter, in accordance with mass-to- light ratios ob- 
served in early-type spiral galaxies ( Moriondo et alJl998h . We 
fix the mass-to-light ratio to 3.80 for the remainder of the fit- 
ting. 

Upon comparing the model rotation curve to the observed 
one we have to decide on which aperture on the two curves 
should have the same radial velocity. Ideally, this corresponds 
to the point which has zero velocity in the rest frame of the 
galaxy, i.e. the point having radial velocity equal to the sys- 
temic velocity. However, because of the high spectral resolu- 
tion of the VISIR data, the systemic velocity of NGC 7582 
known from the literature is not accurate enough. Instead, we 
use the perpendicular spectrum to help constrain the radial ve- 
locity of the dynamical center. If the slit for the perpendicular 
spectrum is aligned along the minor axis of the disk, we ex- 
pect to see a constant velocity along the slit corresponding to 
the systemic velocity. Examining the upper right-hand plot in 
Fig. |3] we see that there is a region over pixels 26-38 with 
constant velocity a;1644 kms -1 . If this corresponds to the sys- 
temic velocity of the galaxy, then the aperture centered on pixel 
40 in the parallel spectrum is the correct center to choose for 
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Table 3. The confidence intervals for the parameters of the 
best-fitting model having M BH - 5.5 x 10 7 M and b =073. 





M BH [10' M Q ] 


b [arcsec] 


68% CI 


[3.2,8.7] 


[0.282,0.341] 


95% CI 


[3.6,8.1] 


[0.273,0.335] 


99% CI 


[4.4,7.1] 


[0.294,0.314] 



the fitting because this is the aperture closest to 1644 kms . 
However, if the region with velocity w 1644 kms -1 in the per- 
pendicular spectrum is affected by an outflow or some bulk mo- 
tion, the center could have a lower radial velocity. The aperture 
centered on pixel 40 in the perpendicular spectrum has a ve- 
locity of 1632 kms 1 , and could therefore be closer to the dy- 
namical center. Comparing again with the parallel spectrum, a 
velocity of 1632 kms -1 lies between the apertures centered on 
pixels 38 and 40. So it is likely that the dynamical center in the 
parallel spectrum either coincides with one, or lies somewhere 
between these two apertures. The most probable systemic ve- 
locity of the galaxy therefore seems to be 1630-1650 kms -1 . 
Correcting to heliocentric velocity (the correction factor at the 
time of observation is -16 kms -1 ), we therefore find that the 
systemic velocity of the galaxy is 1614-1634 kms -1 , approxi- 
mately 60 km s -1 higher than that reported in the Leda database 
based on optical emission lines. Outflows and bulk motion of 
the gas make it difficult to pin-point accurately the aperture cor- 
responding to the systemic velocity. Integral field spectroscopy 
in the near-infrared might help in this respect, although even 
in this case, there may be problems related to dust obscuration 
and outflows. 

By requiring that the observed and the modeled rotation 
curve have the same velocity in the aperture centered on pixel 
38, we derive a best-fitting model with M BH = 1 -Ox 10 8 M Q and 
b =0735 with a reduced^ 2 of 3.76. By shifting the common 
aperture to the one centered on pixel 40, the fit is marginally 
better with a reduced x 1 of 3.36, M B h = 5.5 x 10 7 M and 
b =073. The slightly different best-fit parameters in each case 
serves to illustrate the level of uncertainty related to the diffi- 
culties with knowing the exact systemic velocity. The reason 
we obtain a reduced^ 2 greater than unity is most likely due to 
outflows and bulk motions in the gas that are not included in the 
model. The current modeling is probably the best that can be 
done with these data, and the limitations are due to uncertain- 
ties related to the spatial resolution being larger than the sphere 
of influence of the black hole. 

We choose the model with Mbh = 5.5 x 10 7 M as our best 
fit since this has the lowest^ 2 . This model is plotted together 
with the observed rotation curve in Fig. [5] and a contour plot 
of x 1 is shown in Fig. |6j with 68, 95% and 99% confidence 
regions indicated. The confidence regions on Mbh and b are 
also listed in Table[3] 

5.5. The dynamical center 

The best-fitting value b = 073 indicates that the dynamical cen- 
ter may not have been exactly centered on the slit when the 



1800 r 




1400 r i i i i , 
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Spatial axis [pix] 



Fig. 5. The open circles show the observed rotation curve and 
the solid line indicated the best-fit model with Mbh = 5.5 x 10 7 
M and b =073. The dashed line corresponds to a model with 
b = 0.3 and M BH = 0. 




2x10 7 4X10 7 6x10 7 8x10 7 
log M BH [M ] 



Fig. 6. Contours of^ 2 . The minimum is marked by an 'x', and 
the solid lines mark the 68, 95 and 99% confidence regions. 

spectrum was taken, but instead located 07085 from the edge 
of the slit, west of the center. We found above that the sys- 
temic velocity of the galaxy probably corresponds to the ve- 
locity probed by the apertures centered on pixel 38/40 in the 
parallel spectrum. These constraints can be used to find the lo- 
cation of the dynamical center on the narrow-band image, as 
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6. Discussion 



60 



50 



S 40 



30 



20 



b=0.3" 



20 



30 



40 

y s [pix] 



50 



60 



Fig. 7. The two shorter horizontal lines to the left indicate the 
slit width, and the two vertical lines the center of the two aper- 
tures thought to correspond to the systemic velocity of the 
galaxy. The horizontal line across the image indicates the best 
fit dynamical center position b = 073. The image has been ro- 
tated, so East is up and North is to the right. The axes are la- 
beled according to the slit coordinate system. 



shown in Fig.0 In this figure, the image has been rotated so 
that the length of the slit runs horizontally, hence the abscissa 
is labeled y s . The ordinate is labeled x s since it runs along the 
width of the slit. The location of the slit is indicated by the two 
short horizontal lines on the left-hand side of the image, and 
the line corresponding to b =073 is the horizontal line plotted 
across the full width of the image. The two vertical lines cor- 
respond to the center of the two apertures that were discussed 
in the previous section as likely to correspond to the systemic 
velocity of the galaxy (between aperture 38 or 40). The dynam- 
ical center of the galaxy should therefore lie in the region where 
the horizontal line intersects the two vertical lines. As seen in 
Fig- El this region does not correspond to the exact location of 
the peak in the surface brightness, thereby suggesting that the 
dynamical center may be obscured. 

According to unified models jAntonuccil [T993 ). the loca- 
tion of the dynamical center is expected to coincide with the 
vertex of the ionization cone because the cone morphology is a 
result of the illumination of gas by the AGN continuum, which 
can only escape within + ~ 45° of the axis of the torus. The 
blue outline of the ionization cone is seen in the colour im- 
age in Fig.^ but it is difficult to tell from the image precisely 
where the vertex is. The location of the dynamical center might 
therefore be consistent with the vertex of the cone. 



6.1. Black hole mass estimates in the mid-infrared 

In order to prove the existence of a black hole or other compact 
object at the center of NGC 7582 with spatially resolved gas 
dynamics, the resolution must be sufficient to probe gas dynam- 
ics inside the black hole's sphere of influence. The current data 
do not have the required spatial resolution for this, but if we ac- 
cept that black hole activity is responsible for the AGN activity, 
the data presented here can be used to infer constraints on the 
mass of the black hole. We find a best-fitting black hole mass of 
5.5xl0 7 M , with a 95% confidence interval of [3.6,8. 1]xl0 7 
M , consistent with the prediction of w 4 x 10 7 M from the 
local Mbh-o"» relation. 

This is the first time that a black hole mass in a galaxy 
except our own Milky Way system has been estimated from 
gas dynamics using line emission in the mid-infrared. Usually 
Ha and/or [N n] in the optical have been used, in particular for 
HST data. Ground-based work has so far been limited to, and 
has relied on, emission lines in the near-infrared since better 
resolution is achieved at longer wavelengths. More measure- 
ments are expected to appear in the near-infrared with the appli- 
cation of adaptive optics (Haring-Neumav er et al . 2006). The 
VISIR data presented here are diffraction limited and therefore 
have the best obtainable spatial resolution possible with an 8- 
m telescope at the position of the [Ne n] line. The black hole 
in NGC 7582 seems to be on the limit of what is possible to 
constrain with a spatial resolution of 074, and probably black 
hole masses of ^ 5-10 x 10 7 M in nearby galaxies with gas 
disks can be estimated from the ground using spectroscopy of 
the [Nen]12.8yum line. 

Since dust is easily penetrated at mid-infrared wavelengths, 
gas dynamics in the mid-infrared may prove to be a valuable 
tool for probing the Mbh-o~» relation in galaxies where the nu- 
clear regions are obscured by dust and where optical and/or 
near-infrared spectroscopy may fall short. Gas dynamics in 
the mid-infrared may therefore help determine the black hole 
masses of e.g. Seyfert 2 galaxies where the nucleus is hidden 
behind an optically thick, dusty torus. The assumption is, as 
with optical and near-infrared observations, that the gas is well- 
behaved and has a relatively smooth rotation curve. 



6.2. The morphology of the circum-nuclear disk 

The two bright regions in the disk to the South of the nucleus 
may be dense regi ons of intense s tarformation, deeply embed- 
ded in dust ( Wold & Galliano 2006). The two knots are equally 
bright in [Ne n], bu t one of them is less bright, or undetected, in 
H? an d at 1 1 .9/jm JSosa-Brito et alJEoOlUSiebenmorgen et alJ 
2004). Our modeling shows that the rotation curve of the disk, 
including the two bright knots to the South of the nucleus, is 
consistent with a thin disk rotating in the gravitational poten- 
tial of the bulge. This suggests that the two knots lie in the 
plane of the disk. According to the unified model the ionization 
cone axis is perpendicular to the plane of the torus, and since 
the circum-nuclear disk is probably an extension of the torus, 
the two knots must be shielded from the ionizing continuum of 
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the AGN. We therefore conclude that the AGN continuum can- 
not be responsible for the destruction of PAHs and H2 in one 
of the knots. The lack of PAHs and H2 in one knot relative to 
the other may be explained by differences in the ionizing con- 
tinuum and/or by differences in the physical conditions in the 
molecular clouds. For instance, the PAH- and H2 -deficient knot 
may be further along in its lifetime and has ionized more of its 
local interstellar medium, destroying both PAHs and H2 in the 
process. Alternatively, the PAH-deficient knot may contain a 
larger population of hot O-stars than the other since PAHs are 
found to be better tracers of B-st ar populations rather than very 
hot O-stars dPeeters et alJl2004 . Finally, as PAHs and H 2 can 
arise in different regions and H2 is very temperature-sensitive, 
differences in the physical conditions of the two knots may also 
influence their appearances. 

7. Conclusions 

We have presented data taken with the ESO VLT mid-infrared 
instrument, VISIR, consisting of a narrow -band image centered 
on the [Ne 11] 12.8//m line, and two high-resolution (R = 16000) 
spectra also centered on the [Nen] line, of the composite 
starburst-Seyfert 2 galaxy NGC 7582. The galaxy has a circum- 
nuclear star forming disk and the two long-slit spectra were 
obtained with the slit parallel and perpendicular to the major 
axis of the disk. The parallel spectrum shows a clear rotation 
curve, with an amplitude of * 190 kms _I , and the perpendic- 
ular spectrum, a roughly constant velocity profile as expected 
for a geometrically thin disk. 

We model the disk as a thin circular disk rotating in the 
combined gravitational potential of the bulge and a central 
black hole, and find a best-fit black hole mass of 5.5x1 7 M 
with a 95% confidence interval of [3.6,8. 1]xl0 7 M Q . In the 
modeling we also allow for the fact that the dynamical cen- 
ter may not have been exactly at the center of the slit during 
observation. The best-fitting model indicates that it probably 
was located close to the edge of the slit, approximately 073 to 
the west of the slit center. Taking advantage of the high spectral 
resolution, we argue that the (heliocentric) systemic velocity of 
the galaxy is cz ~ 1614-1634 km s -1 . 

There are two dense starforming knots in the disk, equally 
bright in [Nen], but one shows weaker emission in H2 and at 
11.9/im, indicating that the PAHs and the H2 may have been 
destroyed. Since the rotation in the outer parts of the disk, in- 
cluding the two knots, is well explained by a thin disk rotating 
in the gravitational potential of the bulge, we conclude that the 
two knots lie in the disk and are thereby shielded from the AGN 
continuum. If an ionizing continuum is responsible for the de- 
struction of PAHs and H2 in one of the knots, it must come 
from the star forming region itself rather than the AGN. It can- 
not be excluded however, that different physical conditions, e.g. 
temperature, within the two knots may be responsible for their 
different appearances in H2 and at 1 1.9/um. 

The data are diffraction limited with a spatial resolution of 
» 074, and we argue that this resolution is sufficient to put 
constraints on black holes masses of ^ 5-10 x 10 7 M . Now 
that STIS on board HST is no longer operational, other tech- 
niques for estimating black hole masses in galaxies need to 



be explored. We have shown that gas dynamical estimates of 
black hole masses based on high-resolution mid-infrared spec- 
troscopy are competitive with other techniques, particularly if 
the galaxy nucleus is obscured by dust. Such observations also 
allow us to study obscured star formation in the centers of 
nearby AGN, thus advancing our understanding of the connec- 
tion between AGN activity and nuclear star formation. 
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